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We discuss the evolution of purity in mixed quantum/classical approaches to electronic nonadiabatic dynamics 
CNj ' in the context of the Ehrenfest model. As it is impossible to exactly determine initial conditions for a realistic 

' system, we choose to work in the statistical Ehrenfest formalism that we introduced in Ref. [J. From it, we 

i develop a new framework to determine exactly the change in the purity of the quantum subsystem along the 

' evolution of a statistical Ehrenfest system. In a simple case, we verify how and to which extent Ehrenfest 

ly-^ , statistical dynamics makes a system with more than one classical trajectory and an initial quantum pure 

state become a quantum mixed one. We prove this numerically showing how the evolution of purity depends 

i 1 on time, on the dimension of the quantum state space D, and on the number of classical trajectories N of the 

initial distribution. The results in this work open new perspectives for studying decoherence with Ehrenfest 
ll | dynamics. 



2 of Ref. i). 

In J. C. Tully's surface hopping methods^ for exam- 
ple, the deconstruction goes to BOMD and the recon- 
struction proceeds by allowing the system to perform 
certain specially designed stochastic jumps between adi- 
abatic states. 

In the decay of mixing formalism by D. G. Truhlar and 
coworkers^ the deconstruction stops at the ED and the 
reconstruction is developed by adding decoherence to it. 
This has been shown to be more accurate than surface 
hopping methods for non-Born-Oppenheimer collisions. 

When considered for a single system, ED is a fully co- 
herent semiclassical method, and hence purity preserv- 
ing. As decoherence must be a property of any realistic 
modelj^r— many MQCD models have been reconstructed 
in order to produce electronic decoherence. See, for ex- 
ample, Refs.y-Qjj; that range from one of the most classic 
in this matter— to one of the most recent J^. In Ref. IbH . 
we can find a recent study of decoherence in the context 
of surface hopping and an important conclusion: averag- 
ing over a swarm of initial conditions, decoherence can 
be measured; but the method cannot capture all the ob- 
served effects (for example, the averaging is not enough to 
capture the physics of wave packet bifurcations on multi- 
ple surfaces). In the case of decay of mixing formalisms, 
where the trajectories in the swarm are considered as 
independent, decoherence phenomena are incorporated 
algorithmically (see Eq. 18 of Ref. H). 

In this paper we study the problem from a different 
perspective. First of all, we consider a complete statis- 
tical description of ED, which was introduced in Ref. 1 



I. INTRODUCTION 
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c/2 . The Schrodinger equation for a combined system of 
electrons and nuclei enables us to predict most of the 
chemistry and molecular physics that surrounds us, in- 
cluding biophysical processes of great complexity. Unfor- 
tunately, this task is not possible in general, and approxi- 
mations need to be made; one of the most important and 
successful being the classical approximation for a num- 

• ber of the particles. Mixed quantum-classical dynamical 
' (MQCD) models are therefore necessary and widely used. 

qq ■ We could say that, typically, the technique used to 
build MQCD models is a partial 'deconstruction' of the 

• quantum mechanics (QM) of the total system (electrons 
' and nuclei) followed by a 'reconstruction' that tries to 
, recover the essential properties of the total Schrodinger 

CnI ' equation lost in the deconstruction process. It is unreal- 
istic to expect the reconstructed theory has the same pre- 
dictive power as the Schrodinger equation, so the recon- 
structed theory will apply with enough accuracy only to a 
subset of systems and questions; a subset whose bound- 
aries are difficult to predict a priori. In the literature, 
there are at least two common levels of deconstruction, 
one further away from the total Schrodinger equation for 
electrons and nuclei, called Born-Oppenheimer molecu- 
lar dynamics (BOMD), where electrons are assumed to 
remain in the ground state for all times, and another one 
closer to it, called Ehrenfest dynamics (ED), where nu- 
clei are still classical (as in BOMD) but the electrons are 
allowed to populate excited states (some misleading no- 
tation used in the literature on ED is clarified in Section 
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in full detail. Based on that construction, we develop 
a description of Quantum Statistical Mechanics which 
can be adapted to mixed quantum-classical systems in 
a straightforward but rigorous manner. This allows us 
to consider, in a simple way, the evolution of the purity 
of the quantum subsystem of our Ehrenfest model. We 
prove thus that, while a single Ehrenfest system evolves 
preserving the purity of the quantum state, the behav- 
ior changes dramatically when a statistical distribution 
is considered and, in general, it introduces a change in 
the purity of the quantum system along the trajectory. 
Therefore, we can claim that our general statistical ED 
formalism is purity non-preserving; a property which al- 
ways accompanies decoherence phenomena (see for exam- 
ple Sec. 3.5 in Ref . Il5h . This work opens thus a new line 
of research in the direction of taking into account deco- 
herence phenomena in ED. Of course, those ingredients 
still missing for a proper description of decoherence in 
our Ehrenfest statistical formalism could be later added 
in a reconstruction process similar to the ones mentioned 
before, but this time starting from a, presumably better, 
purity non-preserving dynamics 

A full study of the decoherence process is a very com- 
plex task which includes deep quantum theoretic con- 
cepts as the measurement problem and the interpreta- 
tions of QM (see, for example, Refs. UH and [l6| for a 
general discussion, and Ref. [j for the analysis of the de- 
coherence phenomenon on molecular systems). In this 
paper we voluntarily restrict ourselves to a simple prop- 
erty reflecting the decoherence phenomenon. Namely, the 
change in the 'degree of mixture' of the quantum state 
in a MQCD model, as quantified by the purity Trp 2 . As 
mentioned, we shall see in Sections IIVI and [V] that ED 
provides a framework where this change takes place. The 
actual relation with the electronic decoherence in molecu- 
lar systems requires a much more involved analysis which 
will be developed in the future. 

Besides the incomplete description of decoherence, 
usual approaches to ED have also been often criticized 
on the basis that it does not yield the Boltzmann equi- 
librium distribution for the electrons exactlyji£~— The 
lack of this property, which we agree is desirable, is how- 
ever not enough to rule out ED for all applications, as 
we recently arguedj 22 ' 23 

The structure of the paper is as follows: Sections 
from HIl to IIVI introduce the mathematical formalism and 
the relevant definitions, which are then put into prac- 
tice in the numerical example in Sec. [V] Sec. [TT] reviews 
the notion of purity in QM and proves the well known 
fact that ED preserves the purity of the quantum sub- 
system when we consider the evolution of a single trajec- 
tory from perfectly determined initial conditions. Sec. lIIII 
presents a very brief summary of the formulation of geo- 
metric QM (see [l] for a more careful presentation) and it 
provides an analogous formulation of a quantum statis- 
tical system within the same framework. In particular, a 
suitable formulation of the purity of a quantum system 
is introduced. Sec. IIVI presents the main contribution of 



the paper: first, we review the geometrical formulation of 
ED and its associated statistical equilibrium introduced 
in [J Then, we adapt the tools introduced in the previous 
section in order to be able to study the evolution of the 
purity of the quantum subsystem in a suitable way, and 
to show that ED is purity non-preserving. The use of 
the geometrical formalism, as we will see in what follows, 
allows to perform a very direct analysis of the problem. 
In Sec. [Vj we numerically illustrate the change in purity 
produced by ED using a very simple but extremely useful 
example: a statistical system defined by a pure quantum 
state and an ensemble of initial conditions of the classical 
subsystem. Such a system has been used in the literature 
as a natural framework for molecular dynamics (see for 
example Refs. SEH and[H|). We use it as the simplest 
nontrivial Ehrenfest statistical system where we can show 
how the purity of the quantum part of the system evolves 
in time depending on the coupling between the classical 
and quantum systems, the initial momentum of the clas- 
sical particles, the dimension of the quantum state space 
and the number of trajectories considered in the initial 
conditions. In Sec. I VII we present our conclusions and 
our plans for future works. 

II. PURITY 

A. Purity preservation in quantum mechanics 

Given a Hilbert space H, we shall call density states 

to the elements p obtained as convex combinations of 
rank-one projectors p — {pi, . . . ,pk}, each element satis- 
fying 

= pj, pj = pj, Trpj = 1, j = l,...,k, 

with a probability vector, p :— (pi, . . . ,pk) with Y] ^ pj — 
1 and pj > 0, Vj. The expression of a general density 
state is then 

3 

The evaluation of some observable A on this state is given 
by 

(A)=J2p 3 Mp 3 A)^Tt( P A). (2.1) 

3 

The state of the quantum system is said to be pure 
if the density matrix which represents it is a rank-one 
projector, i.e., if the convex combination above contains 
only one term. If this property does not hold, the system 
is said to be in a mixed state, since, from the physical 
point of view, there is a statistical mixture of the different 
pure states represented by the density matrices pj above. 

Being a Hermitian operator, the matrix p can be diag- 
onalized. Its eigenvalues {Ai, . . . , Xk} satisfy 

< Aj ■ < 1, Vj. (2.2) 



3 



If the state is pure, there is one eigenvalue equal to one, 
the rest being zero. Obviously, the rank of p as a projec- 
tor on the Hilbert space H coincides with the number of 
nonvanishing eigenvalues. 

It is an immediate property that a state p is pure if 
and only if 



TV = !• 



(2.3) 



The proof requires only Eq. f|2 . 2|) . 

The description of a quantum system in terms of a den- 
sity matrix uses von Neumann's equation to introduce 
the dynamics. Then we know that, given the Hamilto- 
nian operator H , the evolution of the state p is given 
by 



ihp(t) = [H(t), P (t)], 



(2.4) 



where [•, •] is the usual commutator of operators. 

Using a simple proof which is formally identical to the 
one that we shall present in the next section, one can 
easily show that this dynamics is purity-preserving, i.e., 
d{Trp 2 )/dt = 0. 



B. Purity preservation in non-statistical Ehrenfest dynamics 

The Ehrenfest equation a 22 ! 23 for a system composed of 
a set of M classical particles (typically nuclei; described 
by the phase space variables R := (Pi, . . . , Rm), P '■= 
(Pi , . . . , Pm)) and a set of n quantum particles (typically 
electrons; described by a wavefunction tp, defined on the 
space parameterized by r := (r*i, . . . , r„)) are: 



(2.5) 



,dH e 



Pj (t) = -{i>(t)\-^(R(t))m)), (2-6) 
oRj 



ih-\m) = H e (R(t))\m), 



(2.7) 



where J = 1, . . . , M and the electronic Hamiltonian op- 
erator H e is related to the molecular one H and it is 
defined as follows: 

h b (r) : =-ft 2 y Iv 2 + T ^ y ZjZk _ 



J_ 1 1 y 

47re ° jtk |fj " ?k 1 47160 Ijfij- fj | 
^E^ 2 , (2.8) 



where all sums must be understood as running over the 
whole natural set for each index, Mj is the mass of the 
J-th nucleus in units of the electron mass, and Zj is the 
charge of the J-th nucleus in units of (minus) the electron 
charge. 



At first sight, given the similarity between Eq. (|2.7p 
and the Schrodinger equation for an isolated full- 
quantum system, one might erroneously think that the 
Ehrenfest evolution for the quantum part of the system is 



unitary' 



If this was correct, then it would be trivial 



to prove that ED is purity preserving, but this is not the 
case. As is well known, for a one to one transformation, 
we can define unitarity as the property of preserving the 
scalar product, i.e., given two arbitrary vectors ip and 
(f>, we say that U is unitary if (Uip\U(p) = (<f>\<t>)- One 
can easily see that any reversible transformation U that 
enjoys this property is necessarily linear: 

(U((pi + ifi2)\U(/)) = (cpi + <p-z\(j>) = {yi\4>) + (y? 2 |0> 

= {u Vl \u4>) + (u<p 2 \U4>) = ((Utp! + u<p2)\Uif>). 

As U is reversible, U<p is an arbitrary vector, and there- 
fore, we must have, 

U((pi + ip 2 ) = U(px + U<p 2 - 

But, although the quantum part of the the equations of 
motion in (|2.7[l resembles a typical Schrodinger equation, 
the coupling with the classical part makes the evolution of 
the quantum system nonlinear. Consequently, it cannot 
be a unitary transformation as defined above. 

Despite this non-unitarity, it is very simple to prove 
that, if we consider the evolution of a single trajectory 
(i?, P, ip) of an Ehrenfest system, the quantum part is 
always in a pure state: 

Theorem 1. Let {R,P,ip) be the initial state of an 
Ehrenfest system subject to the dynamics given by eqs. 
(12. 5p - (I2.7P . Then, the quantum part of the system is al- 
ways in a pure state 

Proof. We consider the density matrix p = \tp) (tp \ corre- 
sponding to the quantum part of the Ehrenfest system. 
The evolution of \ip) is given by Eq. (|2.7p . which induces 
a von Neumann-like evolution for the density matrix at 
every time t 

ihp(t) = [H e (R(t)),p(t)l 
being H e the electronic Hamiltonian in (|2.8p . Then, 
d, 



dt 



Tip 2 = 2Tr(pp) = 2Tr([H e ,p}p) 



= 2 (Tr(P e pp) - Ty(pH e p)) = 0, 
for all times t. Hence, if Tip 2 = 1 at t = 0, it will remain 



so. 



□ 



The main goal of the rest of the paper is to prove that, 
when we consider the case of a statistical ensemble of 
Ehrenfest trajectories, this is no longer the case: the evo- 
lution of an ensemble whose quantum part is a pure state 
at t = will become an ensemble in which the quantum 
part is mixed as long as the initial conditions for the nu- 
clei are not perfectly determined. Thus, in such a case, 
the ED produces a purity change, a necessary condition 
for decoherence. 
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III. GEOMETRIC QUANTUM STATISTICAL MECHANICS 
A. Geometric quantum mechanics 

The aim of this section is to provide a very brief sum- 
mary of the mathematical formalism more thoroughly 
introduced in [l| and references therein. 

Classical mechanics can be formulated in several math- 
ematical frameworks each corresponding to a different 
level of abstraction: Newton's equations, the Hamilto- 
nian formalism, the Poisson brackets, etc. Perhaps its 
more abstract and general formulation is geometrical, in 
terms of Poisson manifolds. Similarly, QM can also be 
formulated in different ways, some of which resemble its 
classical counterpart (see Ref. [13 for a classical reference 
in the context of molecular systems, Ref. [l| for a more 
recent one, and Refs. l28M3ll for more mathematical ap- 
proaches). For example, the observables (self-adjoint lin- 
ear operators) are endowed with a Poisson algebra struc- 
ture (based on the commutators) almost equal to the 
one that characterizes the dynamical variables in classi- 
cal mechanics. Moreover, Schrodinger equation can be 
recast into Hamilton's equations form by transforming 
the complex Hilbert space into a real one of double di- 
mension. The observables are also transformed into dy- 
namical functions in this new phase space, in analogy to 
the classical one. Finally, a Poisson bracket formulation 
has also been established for QM, which permits to clas- 
sify both the classical and the quantum dynamics under 
the same heading. 

This variety of formulations does not emerge from aca- 
demic caprice; the succesive abstractions simplify further 
developments of the theory, such as the step from micro- 
scopic dynamics to statistical dynamics: the derivation 
of Liouville's equation (or von Neumann's equation in 
the quantum case), at the heart of statistical dynamics, 
is based on the properties of the Poisson algebra.- 

Consider a basis {l^fc)} f° r the Hilbert space H. Each 
state \tp) £ H can be written in that basis with complex 
components (or coordinates, in more differential geomet- 
ric terms) {zk}- 

k 

Now, we can just take the original vector space inher- 
ent to the Hilbert space, and turn it into a real vector 
space (denoted as Mq), by splitting each coordinate into 
its real and imaginary parts: 



U 3 z k = q k + ip k h-> (q k ,p k ) e 



n2ra 



Mr 



We will use real coordinates (q k ,Pk), k = 1, ...,n, to 
represent the points of H when thought of as real vector 
space elements. From this point of view, the similarities 
between the quantum dynamics and the classical one will 
be more evident. It is important to notice, though, that 
despite the formal similarities these coordinates (q k ,Pk) 
do not represent physical positions and momenta of any 



actual system. They simply correspond to the real and 
imaginary parts of the complex coordinates used for the 
Hilbert space vectors in a given basis. 

The scalar product of the Hilbert space is encoded in 
three tensors defined on the real vector space Mq. The 
interested reader is addressed to Ref. [l] for the details. We 
just highlight here that two of these tensors correspond 
to a metric tensor g and a symplectic one u) which allow 
us to write the expression of the Schrodinger equation 
as a Hamilton equation, in a form which is completely 
analogous to the Hamiltonian formulation of classical me- 
chanics. It is precisely this similarity the key ingredient 
to successfully combine classical and quantum mechan- 
ics in a well-defined framework to describe the Ehrenfest 
equations (|2.5p - (|2.7p as a Hamiltonian system, as we will 
summarize later and it can also be seen in Ref. [J. 

In this formalism, instead of considering the observ- 
ables as linear operators (plus the usual requirements, 
self-adjointness, boundedness, etc.) on the Hilbert space 
Ji, we shall represent them as functions defined on the 
real space Mq. The reason for that is to resemble, as 
much as possible, the classical mechanical approach. But 
we cannot forget the linearity of the operators, and thus 
the functions must be chosen in a very particular way. 
The usual choice is inspired in Ehrenfest 's description of 
quantum mechanical systems and defines, associated to 
any operator A on "H, a function of the form: 



f A ty) := \{^,A^) 



(3.1) 



The operations which are defined on the set of opera- 
tors can also be translated into this new language. Thus, 
the associative product of operators (the matrix product 
when considered in a finite dimensional Hilbert space), 
the commutator (which encodes the dynamics) and the 
anticommutator can be written in terms of the functions 
of the type defined in Eq. (|3.1[) . As an example, we can 
write the case of the commutator, which will be used 
later: Given two operators A and B, with the corre- 
sponding functions Ja and fg, the function associated 
to the commutator i[A, B] = i(AB — BA) (the imaginary 
unit is used to preserve hermiticity) is written as 



fi[A,B] — {fA, /b} = 



2^\dq k dp k dp k dq k 



> 2 ) 

Thus, from the formal point of view, the operation is 
completely analogous to the Poisson bracket used in clas- 
sical mechanics. 

Another important property in the set of operators of 
QM is the corresponding spectral theory. In any quantum 
system, it is of the utmost importance to be able to find 
eigenvalues and eigenvectors. We can summarize these 
properties in the following result: If Ja is the function 
associated to the observable A, then, as a consequence of 
Ritz's theorem, 32 

• the eigenvectors of the operator A coincide with the 
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critical points of the function f a-, i.e., 
dfA(4>) = ^ is an eigenvector of A. 

• the eigenvalue of A at the eigenvector if) is the value 
that the function j 'a takes at the critical point ip. 

As usual, the dynamics can be implemented in essen- 
tially two different forms (but always in a way which is 
compatible with the geometric structures introduced so 
far): the so-called Schrodinger and Heisenberg pictures 1 . 
In the Heisenberg picture, which is the one we will use in 
what follows, the dynamics is introduced by translating 
the well-known Heisenberg equation into the language of 
functions: 



ifrfA = {/aJh}, 



(3.3) 



being fn the function associated to the Hamiltonian op- 
erator and A any observable. 



B. Geometric quantum statistical mechanics 

1. The probability density and the density matrix 

A classical result in QM states that, given a quantum 
system, the average value of any observable A can always 
be computed as the trace of the observable and some 
density state p, as defined in Section III Al 



(A) = Tr(pA). 



(3.4) 



This result is known as Gleason theorem (see Ref. 33 for 
details). 

Instead of using the density matrix p, we can use an 
alternative approach which is formally closer to the de- 
scription of classical statistical systems and is used, for 
instance, in Ref. [34l Consider a probability distribution 
Fq on Mq and a volume element dpQ, satisfying the 
properties: 



f MQ dfi Q (i>)F Q (i>) = 



= 1. 



• Expected values can be computed as 



(A) 



M ^)W)§f 



(3.5) 



for all /a of the form (|3.1[) : A being a Hermitian 
operator. Notice that we have chosen to integrate 
over all the states in Mq and divide by the norm of 
the state, as it is done in the final section of Ref. [l|. 
This is equivalent to integrate over the states of 
norm one as it was also done in the first sections of 
Ref. H. 

The canonical symplectic form of Mq described in 
Ref. [l] provides a natural candidate for the volume form 
since it is also preserved by the quantum evolution (see 
Ref. [l] for the technical details) . 

Some simple examples for the distribution Fq can also 
be provided: 



Example 1. 

we can use 



For the case of the pure state p 



F Q (i>) = 6(ip-ip ) 



_ I^O><-0o| 



(3.6) 



to satisfy the above two equalities. Analogously, a mixed 

state p = EfcPfc '(tfahM ( where J2kPk = 1 and Pk > 0) 
can be represented by 



Fq(*P)=Y.Pk 6 ^~^ 



(3.7) 



In particular, it is straighforward to prove that, in this 
case, 

Jm q y (w> 

f° r p k = and p = T. k Pkpk- 

The definition of the function Fq contains a number of 
ambiguities which are explained in detail in Ref. [l|. Es- 
sentially, we can add to any Fq a term which integrates 
to zero and has vanishing second-order momenta. Due to 
the structure of the observable functions , this modifi- 
cation will not change any average value computed as in 
Eq. (|3.5p . nor will it change the normalization condition 
for Fq . This defines an equivalence class of distributions 
that produce the same average values, and (through the 
relationship between distributions and density matrices) 
Gleason theorem implies that there is always a distri- 
bution in the class in which the fact that we are dealing 
with a convex combination of rank-1 projectors is visible. 
This is what we used in the example above. 

In order to advance in the formulation and illustrate 
these facts more precisely, we can consider, for every 

) G H, the following function: 



(3.8) 



Now, it is easy to see that the following (averaged, now 
i/i-independent) function 



f P (v) 



(3.9) 



is the phase-space function associated to a density oper- 
ator p defined by 



f P (v) = ■xivlPV)- 



(3.10) 



(3.11) 



With these definitions, we can now prove the following 
result: 
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Theorem 2. Let us consider a quantum system in a 
state described by a probability distribution Fq . Then, the 
equivalence class of distributions that produce the same 
expected values as Fq contains as an element: 



(3.12) 



where {ipk} *s the set of critical points of f p (as defined 
in Eqs. iTO|) and iTOT]) ) and 

Afc = fp(ipk)- 
Proof. It is immediate if we realize that 

(A) = I d» Q &)F Q &)^l 
Jm q (w> 




Tr(pA). 
(3.13) 



Indeed, the operator p appearing in this expression and 
defined in Eq. (|3.10[) can be shown to be a density ma- 
trix (i.e.. p 2 = p. p + = p, and Trp = 1). and Gleason's 
theorem guarantees it is unique. 

We also know that, if we use the spectral decomposi- 
tion of p, i.e., p = J2k X kPk, with p k = ^jgfffi , being ip k 
and Afc its eigenvectors and eigenvalues, respectively, we 
also have that 



= ^A fc Tr( Pfe A) =Tr(pA), 



as we set out to prove. 



fA^k) 

(4>k\4>k) 



(3.14) 



□ 



Hence, from Gleason theorem, we know that, among 
all the equivalent distributions, there is always one Fq 
equal to a convex combination of Dirac-delta functions. 
Notice that the function Fq provides us with all the infor- 
mation encoded in the density matrix p. As a probability 
density, it allows us to define the average values of the 
observables, and in the form f p , it allows us to read the 
spectrum of p from the set of critical points. 

This result allows us to realize in terms of Fq any quan- 
tum system: as the average values coincide with those 
obtained from the spectral decomposition of the density 
matrix, we can use it to implement any desired model. 
We will see a practical example in the following sections. 



2. Geometrical computation of purity 

Finally, we would like to analyze purity in this geomet- 
rical context. We saw in Section [FT] that purity preser- 
vation is encoded in the behavior of p as a projector. 
If the evolution of the system preserves the purity of 



the density matrix we say that the evolution is purity- 
preserving, while in the other case, we call it purity 
non-preserving. 

Our first task is to express in this geometrical language 
the concept of purity. Consider the following expression 
for any f p (ip): 



(P) - 



(3.15) 



Then, (p) — 1 if the state is pure and (p) < 1 if the state 
is mixed. This is so because, trivially, 

(p) = Tr(p.p) = Trp 2 . 

Let us now mention how this change in the purity of 
the state can be detected in the measurement of aver- 
age values of observables. Recall that a change in the 
purity as above produces a transformation at the level 
of the states such that a pure state corresponding to a 
distribution of the form Fq(i/j) = S(ip — tpo) becomes a 
distribution of the form Fq(t/)) = PkS(ip ~ V'fc)) where 
there is more than one value of pj different from zero. 
Then, it is immediate to prove that the average value of 
a generic observable A will be different between one case 
and the other. 



IV. PURITY CHANGE IN EHRENFEST STATISTICS 
A. The definitions 

In this section, we will now extend the previous con- 
struction to the Ehrenfest case, by combining it with the 
approach introduced in Ref. [l|. 

First, let the physical states of our Ehrenfest system 
correspond to the points in the Cartesian product 



M = M c x M Q , 

where Mq is the phase space of the classical system. The 
physical observables will be now functions defined on that 
manifold. To define statistical averages of observables 
depending on classical and quantum degrees of freedom 
(i.e., functions as /a^jVO) we consider 



M c xM Q (V#) 



(4.1) 



where £ = (R, P) £ Mq represents the classical degrees 
of freedom, ip = ip(q,p) € Mq the quantum ones, and 
d/J'QC — dpQdpc is the volume on the state space mani- 
fold M. 

We can now ask what properties we must require from 
Fqc in order for Eq. (|4.ip to correctly define the statis- 
tical mechanics for the Ehrenfest dynamics (ED). Anal- 
ogously to what happened in the quantum case, the con- 
ditions are as follows: 
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The expected value of any constant observable 
should be equal to that constant, which implies that 
the integral on the whole set of states is equal to 
one: 



M c xM Q 



dti Q c{i,ip)F QC {i,iP) = l. 



(4.2) 



• The average, for any purely quantum observable Ja 
of the form in (I3.1[) , associated to a positive definite 
Hcrmitian operator A, should be positive. This im- 
plies the usual requirement of positive probability 
density in standard classical statistical mechanics. 

In Ref. [jj, it was proved that the ED defined on the 
manifold M is Hamiltonian with respect to the Poisson 
bracket 



{•> -}qc 



{•) -}c + '}q> 



(4.3) 



where {•, -}c represents the usual Poisson bracket of the 
classical degrees of freedom, and {-,-}q represents the 
Poisson bracket defined by Eq. (|3.2[) . 

Being Hamiltonian, we know that we can define an 
invariant measure on the space of states M. We shall 
denote such a measure by dfj,Qc- Thus, the dynamics 
defined on the microstates is straightforwardly translated 
into the probability density Fqc as a Liouville equation: 



Fqc = {Ih, Fqc}qc, 



(4.4) 



where fu is the Hamiltonian function of the Ehrenfest 
system: 

fmp , \- P/ , (i>(q,p)\He(R)mq,p)) 
P; ^■=T^2M J + W? ,p Mg ,p)> ■ 

(4-5) 

Analogously, the evolution of any function /a(£>V0 is 
given by 



Again, this property provides us with a natural can- 
didate for the volume element dfj,Qc (and, for analogous 
reasons, also for dpLc) arising from the symplectic form 
which gives its Hamiltonian structure to the Liouville 
equation in this context. As it happens in the pure quan- 
tum case .this volume form is preserved by the dynamics 
(see Ref. 13). 

With this in mind, we can consider the analogue of the 
objects introduced in the previous section. Hence, given 
an Ehrenfest system in a state described by a probability 
density Fqc(£,T]), where £ = (R,P) and r/ — r)(q,p), we 
can consider the definition of the operator 



(4.7) 



which still depends on the classical variables £ and, there- 
fore, it can be interpreted as having turned the phase- 
space representation of the quantum part into the more 



familiar one based on density matrices. Also, since we 
have not integrated over £, this object still represents in 
a certain way a probability density in the classical part 
of the space. 

As for any £- dependent operator (see Ref. [l|), we can 
define the associated phase-space function: 

m,V):= J^^QC{^)^^. (4.8) 

It is also possible to integrate again the object in 
Eq. l[4T7) . and define: 



P 



d»(Op(0, 



(4-9) 



M c 



which is a purely quantum object encoding the averaged 
information of the complete system. Notice that, as it 
is usually done and in order to lighten the notation, we 
will often use the same symbol for different objects (as 
in p(£) and p), understanding that it is the explicit indi- 
cation of the variables on which they depend what dis- 
tinguishes them notationally. Also, for simplicity, we use 
the same symbols for operators in the full-quantum case 
in sec. IIIIB1 and for the ones in the quantum-classical 
scheme in this section. 

We can now formulate the dynamics in terms of this 
operator. After a brief computation, we obtain: 



dt P = lH 



dn Q c{t,iP)F QC {^) 

MqxM c 



fA(£,1>)={fA(£,1>),fB{Z,ll>)}QC- (4-6) tem\ij ),i.e 



(V#> 

(4.10) 

where H e (£) is the very electronic Hamiltonian defined 
at the beginning of the paper. In ED, considered statisti- 
cally, this equation represents the analogue for the mixed 
case of von Neumann's equation. 

Example 2. If we consider a single ('pure') state of the 
classical system £o an d a pure state of the quantum sys- 



F Q c(M = S^-^)S(i>-il>o), 



we obtain 



\(il>o,v)\' 

(-00 I ^o) 



Analogously, the associated operator in Eq. (|4.7p reads 

l^oX^oi 



and the one in Eq. (|4.9p is 



= l^oMol 

as expected. 

We can also consider the two marginal distributions: 
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The distribution in Mq obtained by integrating out 
the classical degrees of freedom: 



F Q ty) = / diic(OF QC (£,il>), (4.11) 
Jm c 

and the corresponding classical version in Mc- 



Fc(0= / dp Q ^)F QC {^^). (4.12) 
Jm q 

Obviously both functions are distribution functions on 
the corresponding manifolds with analogous properties to 
Fqc , and they could be used to compute expected values 
of functions depending only on ip or £, respectively. 

Example 3. It is immediate to check that the definitions 
make sense for distributions of the form 



(4.13) 



i.e., for 'pure' classical part and a quantum-mixed one 
canonically expressed with deltas. 

In this case, the marginal distributions are of the form 

F Q (VO=$>fc<KV>-Vfc), Fc(0 = *(£-£o). (4-14) 



Also notice that, in terms of the quantum marginal dis- 
tribution Fq , we can write the density matrix in Eq. (|4.9|) 

as 



(4.15) 



Once we have recovered the needed ingredients, we can 
discuss the quantum purity of a system governed by ED: 

Definition 1. We say that a quantum-classical system 
is quantum-pure if and only if 



TV = {p)= f dflQ^FgOj) 
JM Q 



'(V#> 



= 1, (4.16) 



being p the one defined in eqs. (14. 9p and (|4.15p . 

In case that the state of a system does not satisfy the 
condition above, we say that it is quantum- mixed. 

For the sake of completeness, and in order to better 
connect with the purely quantum case we discussed in 
sec. lIIIBl we can consider now the function obtained from 
fp in Eq. (I4.8[) averaging directly over Mc, i.e.. 



f P (ri) 



Mc 

m c xm q wm 

a (mi? t nMV'XV'h) 
&VqW)FqW) — TTTR ■ 

Mq W 



(4.17) 



This function plays the role of the function in (|3.9[) in 
the pure case. Indeed, we have that 

fM = (V\ ( / dp QC F QC ((i, WtTTJy) If) 

\Jm c xMq WW J 

(4.18) 



= fol(7 

\Jm 



and hence it corresponds to the quantum expected value 
of the operator p in eqs. f)4. 9[) and (|4.15p . whose full- 
quantum analogue is the one in (|3.10l) in sec. MI Bl 



B. The application: transferring uncertainty between the 
classical and quantum parts 

Consider the following initial distribution evolving un- 
der ED: 

F Q c(0) = 6((-£o)5(ip-ip ). 

This system is completely deterministic and, therefore, 
the Liouville equation will produce exactly, as a solution, 
the integral curves of ED (£(t), ip(t)). Thus we can write: 

F QC (t) = F QC m, m) = - w - m)- 

Consider now a slightly more complex system, consi- 
tuted by a distribution of N equally probable classical 
states, and a pure quantum state at t = 0: 

Fqc(0) = E S <£ ~ *o)J W - fh). (4-19) 

We can also write the marginal distributions as we did 
in Section IPvTA"! 

i N 

^(0) = T7 E *(£ - ; F Q (0) = 5& - ^o). (4.20) 



k=l 



The evolution of such a system becomes 
1 N 

F Qc(t) = jj E^- $ f&^o;t))*-s;Ko',fo;t)), 



fc=i 



(4.21) 

where ($?(£oS^o;*) ? $^(£oS^o;*)) represents the Ehren- 
fest trajectory having (£oi"0o) as initial condition. 

The evolved marginal distribution in the quantum 
manifold is now: 

1 N 

F Q® = ^E^-^o^o;*)), (4.22) 



k=l 



and then, using Eq. (|4. 16[) . we have: 



A' 



N E /<£* (tk 



(4.23) 
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where, of course, 

/p(0) = |(^o))| 2 . (4.24) 
The associated density matrix at time t reads: 

N h <*;(^^!*)i*;(^^;*)>' ' 

and at time t = 0: 



P(0) - (Vol^o) 
Now, the purity at time t = is 



(4.26) 



(p(0))= / d/x g i?o(0)/ p (0) = 1, (4.27) 



but at a general time f, 



(P(*)> 



dfi Q F Q (t)f p (t) 



N 

E 



(*J(&to;t)l*$(&ito;*)> s 



l 2 ll*S(^,^;*)ll a 



(4.28) 



Hence, the purity of the system seems to evolve in time, 
in general in an involved way. We can compute, though, 
its time derivatives, in order to get an idea about its ini- 
tial evolution. After some calculations which we detail 
in appendix [B] one can see that, for the initial state con- 
sidered in this section, we have 



= 0, 



(4.29) 



and 



4=0 



(4.30) 



which is negative definite unless the expectation value of 
if e (£o) — ^e(Co) a ^ IV'o) vanishes. This means that the 
purity evolves from its initial value, and it does so by 
decreasing, which is entirely expected if we realized that 
we started from a state in which the purity is maximal. 

Thus, we can see that the evolution has made a 
quantum-pure system become a quantum-mixed one, and 
we can claim that: 

Theorem 3. The statistical Ehrenfest evolution is purity 
non-preserving for Hamiltonians with quantum- classical 
couplings. 

This is the main result of our paper: even though the 
ED of a single-trajectory state does preserve purity, when 
we consider a statistical state the behavior changes. In 



the case above (and as we will see numerically in the 
next section), we show how it is possible to transfer un- 
certainty from the classical domain into the quantum one. 
Analogously, it is straightforward to see that an analo- 
gous process happens when we consider a single classical 
state coupled to an ensemble of quantum states. This 
uncertainty transfer takes place whenever the dynamics 
of the quantum and the classical subsystems are coupled 
to each other, since such coupling produces a splitting of 
the trajectories from the given initial conditions and thus 
the mixing of the final states. 



V. NUMERICAL EXAMPLE ON A SIMPLE EHRENFEST 
SYSTEM 

To illustrate the concepts introduced in the previous 
section and to see the purity change in a complex nu- 
merical case, let us consider a situation with an equally 
probable initial distribution as in Eq. (|4.19p of classical 
particles with N — 5 and a quantum part constituted by 
a 10-level system. Such a system has been used in the 
literature as a natural framework for molecular dynam- 
ics (see for example Refs. 0, and l2lh . We consider 
a Hamiltonian function for the quantum-classical system 
of the following form: 



h = J + (ip\A + eJ cos 6B\ 



(5.1) 



where (9, J) are canonically conjugated classical variables 
and A and B are Hermitian matrices acting on the quan- 
tum vector space C 10 . The classical part is written in 
action-angle coordinates to simplify the analysis. 

The dynamics of the system is obtained from the solu- 
tions of the N different trajectories with initial conditions 
defined by each one of elements of the initial classical dis- 
tribution {9 1 (0),9 2 (0),e 3 (0),6 4 {0),9 5 (0)}. The resulting 
distribution takes thus the form given by Eq. (|4.21l) . As 
it can be found in the Supplementary Material, for all 
the trajectories presented below the initial conditions for 
the classical subsystem are chosen as 

0i(O) = 0.9766548288669266, 
2 (0) = 0.5013694871260747, 
3 (O) = 0.9052199783160014, 
4 (O) = 0.5068075140327187, 
5 (O) = 0.9543157645144570 . 

The points were fixed as five random points in the inter- 
val [0,1]. Choosing different initial conditions leads to 
equivalent results. 

The quantum initial condition is chosen to be, for all 
five trajectories, 

= (1,0,..., 0). 

The evolution defines thus a equiprobable dis- 
tribution of N quantum-classical single trajectories 
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^o; t), ^b; t)) e M c x M Q , For each tra- 

jectory, the equations of motion are of Ehrenfest type 
and they are given by eqs. (|2.5l) - (|2.7p . with H e = A + 
eJ cos OB. 

These dynamical equations exhibit two different 
regimes: 

• For e = or J(0) = 0, the classical and the quan- 
tum subsystems evolve uncoupled. The purity of 
the quantum subsystem is always equal to one. 

• For non-zero coupling constants and classical mo- 
menta, the behavior of the system depends sharply 
on the initial conditions: the evolution from differ- 
ent initial conditions for the classical subsystems 
is very different. The purity tends then to evolve 
to the purity corresponding to a set of five ran- 
dom projectors on 1-dimensional subspaces of the 
Hilbert space. 

The evolution of the purity of the resulting system is 
obtained from Eq. (|4.28p after integrating the dynamics 
numerically. Some interesting observations can be ex- 
tracted from the results: 

• In fig. [1] we represent the evolution of the purity for 
a fixed value of the coupling constant e and increas- 
ing value for the initial condition of the classical 
momentum Jo. We see how this makes the sys- 
tem change its originally integrable behavior and 
become more and more chaotic. 




FIG. 2. Evolution of the purity for Jo = 0.8 and e = (dashed 
green line), e = 0.1 (blue line), and e = 0.2 (red line). Again, 
the reference (black line) represents the level of purity of a 
distribution of N = 5 random projectors on C 10 . 



line), in which its value rapidly decreases to a given low 
one and fluctuates around it after that, can be explained 
as a consequence of the dynamics mentioned before. If 
the system exhibits sensitive dependence on the initial 
conditions, after a small lapse of time, the different tra- 
jectories become decorrelated from one another in the 
Hilbert space. Therefore, the density matrix becomes the 
normalized sum of N rank one projectors chosen at ran- 
dom. In this situation, it can be shown (see appendix |A"]) 
that the purity is distributed around an expected value 




FIG. 1. Evolution of the purity for e = 0.1 and Jo = (dashed 
green line), Jo = 0.3 (red line), Jo = 1 (blue line) and Jo = 1.5 
(brown line). We also depict the reference (black line) of the 
level of purity of a distribution of N = 5 random projectors 



• Instead, we can consider a fixed value of the ini- 
tial classical momentum and increase the value of 
the coupling. It can be remarked that the system 
reaches the level of purity of the set of random pro- 
jectors much faster than in the previous case, (see 
fig.©. 

The interesting behaviour of the purity shown for some 
values of the parameters in figs. Q] (brown line) and [5] (red 



E[Tr(p 2 )] = 



N + D-l 

ND : 



(5.2) 



which is represented by the black straight lines in the 
figures, with fluctuations of size 



(E[Tr(p 



2\2l 



E[Tr(p 2 )] 2 



1/2 



V2 / (JV -!)(£>- 1) 
ND V N(D + 1) 



1/2 



(5.3) 



where N is the total number of trajectories used and D is 
the dimension of the quantum Hilbert space. Notice that 
if D is very large, the expected value for the purity tends 
to its minimal value 1/N and the fluctuations tend to 
zero. Remember that the degree of mixture of a quantum 
system ranges from the pure state case (i.e., purity equals 
to one), for which the density operator is a projector 
on a one dimensional subspace of the Hilbert space and 
the maximal mixture case (thus minimum purity) which 
corresponds to a density operator which is proportional 
to the identity matrix. As it must have trace equal to 
one, the proportionality factor is equal to the inverse of 
the dimension of the Hilbert space. 

Figs. [1] and [3] show also the effect of the coupling be- 
tween the classical and the quantum subsystems, as well 
as the effect of the momentum (or equivalently, the en- 
ergy) of the classical particles on the evolution of purity. 
The greater the strength of the coupling and the energy 
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of the classical particle, the faster the system reaches its 
asymptotic behavior. Although clearly the coupling has 
a stronger effect. Thus, analogously to what happens 
in the case of molecular dynamics, where the velocity 
of the classical nuclei induces the coupling between all 
the eigenstates of the electronic Hamiltonian, we see here 
that it also has the effect of mixing the quantum part if 
the Ehrenfest system is treated statistically. 

From our analysis and example, we can then conclude 
that statistical Ehrenfest dynamics provides a framework 
in which the evolution affects the quantum dynamics in 
at least one way decoherence does. Changing the degree 
of mixture of the quantum state is certainly one of the 
most relevant effects of electronic decoherence on molecu- 
lar systems, although further work is required to analyze 
whether or not other decoherence effects can be explained 
by our construction. 



well as on the coupling and the velocity (or energy) of the 
classical system. 'How much' decoherence the statistical 
ED contains, as related to interesting practical applica- 
tions or to other mixed quantum-classical schemes) 36 ' 37 is 
a complex and important question that we shall explore 
in future works. 
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VI. CONCLUSIONS AND FUTURE WORK 

An appropriate description of the electronic decoher- 
ence in molecular systems is on the wishlist of every 
quantum-classical dynamics scheme. In which amount 
each theoretical model includes the sought effects is a 
complicated question whose answer will depend both on 
the model and on the intended application. Loosely 
speaking, we could expect the different models to range 
from 'no electronic decoherence at all' (e.g., BOMD) to 
'a perfect description of quantum electronic decoherence' 
(say, full quantum dynamics of electrons and nuclei) , with 
most of them lying somewhere in between the two ex- 
tremes. Several methods^— have been developed to deal 
with the problem but, to our knowledge, without a defini- 
tive solution. 

Pure Ehrenfest dynamics (ED for single states) clearly 
lies in the 'no decoherence at all' side of the spectrum, 
since it does not even allow the change of purity. How- 
ever, if we are willing to accept the possibility that the 
electrons evolve into a mixed state, i.e., a statistically 
uncertain electronic state, then, to perform a coherent 
analysis, we should also allow the initial conditions of the 
nuclei to be described statistically. This point of view was 
used in some of the papers mentioned above when con- 
sidered within surface-hopping or decay-of-mixing for- 
malisms. In this work we implement the idea using the 
statistical description of Ehrenfest dynamics introduced 
m Ref . LH combined with a specially convenient geomet- 
ric formalism we introduce. Using these methods, one 
can observe that a system starting from uncertain nu- 
clear initial conditions and a pure electronic state evolves 
into a situation in which the electronic state is no longer 
pure but mixed, i.e., ED is capable of transferring uncer- 
tainty between the nuclei and the electrons. Besides, our 
method provides us also with tools to compute in a simple 
but rigorous way the purity change, and it indicates an 
interesting dependence of the effect on the number of tra- 
jectories N and the dimension of the Hilbert space D, as 



Appendix A: Purity of a sum of random projectors of rank one 

In this appendix we derive the expected value of the 
purity and its fluctuations when the density matrix is 
obtained from the sum of decorrelated, random, rank- 
one projectors. 

The purity for a density matrix of the form 



1 N 
AT 51 



IS 



where we have defined 
Xjk ■= 



(i/ij\i/)j){ipk\^k 



Next, we will determine the probability density for X12 
when ipi and ip2 are two decorrelated random vectors. 

Given the global U (N) symmetry of the problem we 
can take the first vector to be 

tp! = (a, (),••• ,0), ael, 

and the second chosen at random. If we denote 

ip2 = (<7i + ipi, 1i + ip2, ...,qn + ipn), 

we get 

X ■— Xl2 — n ' 



where r 2 := Yljilj +Pj) 
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We also need the adequate probability measure in Mq 
that distributes the random vector ip2- It can be defined 
by 

d^Q = f(r 2 )dp Q , 
where / is a positive function chosen so that 



dlTr 



1. 



As we shall see, the actual form of / is not relevant for 
the distribution of the purity. 

It will be convenient to write the probability measure 
in Mq = R 2 x R 2D ~ 2 in the following way: 

cLitq = f(z 2 + R 2 )zdzd8R 2D - 3 dRdn 2 D-3, 

where (z, 9) represents polar coordinates in the plane 
(<?i,Pi), R is the radial coordinate in M. 2D ~ 2 (i.e., R 2 := 
2j>i(Qj +Pj))i while the volume element dCl^DS stands 
for the angular coordinates in M. 2D ~ 2 . In these coordi- 
nates r = V z 2 + R 2 and x — z 2 jr 2 . 

The next step is to perform the change of variables 
from (i?, z) to (r, x) ■ Taking into account that the Jaco- 
bian is 

r 

J = 



one obtains 



and marginalizing out all variables except x we nn d the 
needed measure: 

dn x = (D-l)(l-x) D - 2 d X . 

Once we have determined the probability distribution 
for x we can compute the average value for the purity. 
Using the expression at the beginning of this appendix, 

E[Tr(p 2 )] = ^_ fiV + 2^E[xdl = ^ (l+(AT-l) E[ X ] 

(Al) 

where we have used that all random variables Xij are 
identically distributed. 
Finally, given that 

f 1 1 

E[x] = J xdn x = — , 



we obtain the sought result: 



E[Tr(p 2 )] 



N + D-l 
ND 



As for the fluctuation, one has 
4 



a 2 = 



N 4 



j<k 



2(^-1)^,2. _. l2 ,_ 2(JV-1) (D-l) 



N 3 



W]-E[xT) 



iV 3 D 2 {D + \y 
(A2) 



where we have used that 
E[X 2 ] 



D(D + 1)' 
Appendix B: Derivatives of the purity 



In sec. IIVB1 we explained how an initial uncertainty 
in the classical part of the state of an Ehrenfest system 
produces a change in the purity at t = 0, even if the ini- 
tial state is quantum-pure. In this appendix, we present 
in more detail the calculations that led us to that con- 
clusion. 

First of all, we need the first time-derivative of the 
purity. We know that the evolution of the p(t) is given 
by Eq. (|4.10p . Then, the evolution of the purity can be 
written as: 

j t {p{t)) =2Tr(p(t)p(i)) 

= J dpdp'FF'Tv^hWP^We}- P r ) 

=ih~ l J diidn'FF"Tc([P^Pv\ ■ H e ), (Bl) 

where the integral is taken over (Mq x Mq) 2 , and we 
have denoted 



dp 




(B2a) 


dp! 


= dp Q c(£',ip'), 


(B2b) 


F 




(B2c) 


F' 


= F QC {i l ^'l 


(B2d) 


H e 


= H e (0, 


(B2e) 


K 


= H e (a 


(B2f) 


p^ 


WW 


(B2g) 


iV 


WW) ' 


(B2h) 



In the last step, we also used that 

Tr([P^, H e ] -P^)= Tr(P 4 , ■ H e ■ P 4 , - H e ■ P. 4 , ■ P r ) 

= Tr(H e ■ i> -Pi, -H e -P^,- PW 

= Tr (He -[P fl PW) 
= Tr([P i ,,P i , / ]-H e ). 

Also, as it is common in statistical dynamics, we can 
assign the time-evolution of the state to the probability 
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distribution Fqc and see the objects £, £' and \ip), \ip') 
just as the initial conditions, or we can alternatively think 
that Fqc is the static distribution of initial conditions 
and consider that the time-evolving objects are £, £' and 
\ip), \ip')- Either dynamical image is valid, and the two 
of them produce, of course, the same result, but we have 
performed the calculation thinking in the second way, 
which looked to us slightly more direct. 

Now, in our example, P^(t = 0) = P^,'(t = 0) = P^ . 
Then, we see that the commutator [P^p,P^i] vanishes. 
Thus, we can conclude that: 



!«*» 



t=o 



(B3) 



Using again Eq. (|2.5p - (|2.7p . we can also compute the 
second derivative of the density matrix, 

p(t) = J d^F (r 1 ([{ff e ,/ ff }c,P*]) 

-h~ 2 ([H e ,[H e ,P^]])), (B4) 



where we have used the same notation in eqs. (jB2[) . and 
the integral is this time extended to Mq x Mc- With 
this expression, we can compute the second derivative of 
the purity: 



^{p(t))=Tr(p(t)p(t) + W)y 



(B5) 



Now, using eqs. (|4.10p . (|Bip and (|B4[) we can calculate 
a more explicit form for this second derivatives in terms 
of the objects associated to the geometric formalism 



dt 2 



(p(t)) = I dfidp.'FF' \ ih^TY[\{H e ,f H },P^\ ■ P 4 ,, 
-rr 2r Ir([H e ,[H e ,P^]] 

- h- 2 Tr([H e ,P^] ■ [H' e ,Pv]) 
dp,dp!FF' \ih- lr Tr([{H e J H },P^\ ■ 
+ ?T 2 Tr([ff e , P,p] ■ [H e , (B6) 

- h- 2 Tr([H e , P 4 ,] ■ [H' e ,P r ]^j 
dfidfi'FF' {ih- lr Tr([{H e J H },P^\ -P^ 
+ /T 2 Tr([ff e , P 4 ,] ■ [{H e - H' e ),P^,] 

where we used that 

Tr([HeAHe,Pi,]} -Pp) =Tr([i7 eiJ R ] ■ [# e , -?>])■ 

Using this expression, it is finally straightforward to 
compute the second derivative at t = for the distribu- 



tion Fqc given by Eq. 1)4.2 ip : 

= -h- 2 ^o\(H e (e n )-HM)) 2 Ho) 



(B7) 
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